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Evidence for the dissipation region in magnetotail reconnection 

12 3 

Seiji Zenitani, Iku Shinohara, and Tsugunobu Nagai 



Signatures of the dissipation region of collisionless mag- 
netic reconnection are investigated by the Geotail spacecraft 
for the 15 May 2003 event. The energy dissipation in the 
rest frame of the electron's bulk flow is considered in an 
approximate form D*, which is validated by a particle- in- 
cell simulation. The dissipation measure is directly evalu- 
ated from the plasma moments, the electric field, and the 
magnetic field. Using D*, a compact dissipation region is 
successfully detected in the vicinity of the possible X-point 
in Geotail data. The dissipation rate is 45pWm~ a . The 
length of the dissipation region is estimated to l-2d' oc (lo- 
cal ion inertial length). The Lorentz work W, the work rate 
by Lorentz force to plasmas, is also introduced. It is posi- 
tive over the reconnection region and it has a peak around 
the pileup region away from the X-point. These new mea- 
sures D* and W provide useful information to understand 
the reconnection structure. 



1. Introduction 

Magnetic reconnection is one of the most important pro- 
cesses in many plasma systems. It drives various explosive 
events such as solar flares, substorms in the Earth's magne- 
tosphere, and disruptions in laboratory devices. It is well 
known that the reconnection process is crucially influenced 
by the dissipation physics near the reconnection point (X- 
point). The structure of the dissipation region as well as 
the accommodated physics is of critical importance for the 
understanding of the reconnection mechanism. Significant 
efforts have been made on this subject by theories, numerical 
simulations, and in situ satellite observations. 

In a collisionless kinetic plasma, typically at a distance 
of the ion inertial length from the X-point, ions decou- 
ple from magnetic field lines while electrons remain mag- 
netized. As a consequence, Hall physics plays a role inside 
the ion-decoupling region. Signatures of Hall physics such as 
quadruple magnetic field, bipolar normal electric field, and 
the relevant current loops [Sonnerup, 1979] were confirmed 
by in situ observations in the terrestrial magnetosphere [Na- 
gai et al., 2001; Oieroset et al., 2001]. 

Deep inside the ion-decoupling region, there is a thin layer 
where electrons depart from field lines. A compact dissipa- 
tion region is located in a close vicinity of the X-point. In 
addition, recent simulation works suggest that narrow elec- 
tron jets extend from there, stretching the nonideal layer in 
the outflow directions [Karimabadi et al., 2007; Shay et al., 
2007] . These signatures were recently confirmed by satellite 
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observations. Nagai et al. [2011] found a compact region 
with an intense cross-tail current and neighboring electron 
jets at a super- Alfvenic speed in a magnetotail reconnection 
event. 

Until recently it was not clear how to identify the dissi- 
pation region inside the electron nonideal layer. Taking the 
energy transfer into account, Zenitani et al. [2011a, b] have 
proposed a general measure of the dissipation region. Us- 
ing particle-in-cell (PIC) simulations, they distinguished the 
dissipation region from the outer electron jet, which turned 
out to be an oblique projection of a non-dissipative current 
sheet [Hesse et al., 2008; Klimas et al., 2012]. 

In this Letter, we show in situ observational evidence for 
the dissipation region in a magnetotail reconnection event. 

2. Observation and Analysis 

On 15 May 2003, Geotail detected a reconnection event 
in the premidnight tail. This event was extensively analyzed 
by Nagai et al. [2011] and so we only repeat important fea- 
tures here. Figure 1 shows key properties from 1053:00 UT 
to 1058:30 UT on 15 May 2003. All data are presented in 
the spacecraft (SC) coordinates, because the electric field 
data are available only in the SC coordinates. During the 
period of our interest, the SC coordinates are virtually the 
same as the GSE and GSM coordinates. The rotation an- 
gle to the GSE coordinates is less than 2.4° and the an- 
gle to the GSM is less than 5°. The satellite position was 
(X gsm ,Vgsm,^gsm) = (-27.8, +3.3, +3.57? E ) at 1056 UT. 

Figure la shows magnetic field B z at 16 Hz, measured 
from the magnetic field experiment (MGF) [Kokubun et 
al., 1994]. It suddenly turned from southward (B z < 0) 
to northward (B z > 0) at 1055:44 UT, as indicated by the 
red dotted lines. Judging from this and many other signa- 
tures [Nagai et al., 2011], Geotail encountered a potential 
X-point (hereafter X-point) at this time. 

Other quantities are shown at 12 s time resolution. Fig- 
ure lb shows electric field data, obtained from the electric 
field detector (EFD) [Tsuruda et al., 1994]. Technically, 
they were measured every 3 s and then averaged over 12 
s. We note that the average electric field during 1057:45 - 
1057:58 UT (the arrow in Fig. lb) is represented by one 
during 1057:45 - 1057:52 UT, because EFD had violent, 
high-frequency noises of AE ~ 50mVm _1 during 1057:52 
- 1057:58 UT. One can see that E x has a bipolar signature 
across the X-point. This is a well-known signature of kinetic 
reconnection (e.g., Hoshino [2005]). 

Figures lc and Id show the ion and electron bulk ve- 
locities. Plasma moments are calculated from distribution 
functions that are measured by the low-energy particle ex- 
periment (LEP) [Mukai et al., 1994]. Due to the low time- 
resolution of 12 s, these moments contain field-aligned flows 
along the separatrices from 1055:30 UT to 1057:00 UT. To 
rule out the field-aligned Hall currents, transverse veloci- 
ties are presented (vij,vi e ). The v± ey is distinctly fast 
near the X-point (Fig. Id). As seen in Figure lc, both 
v± ix and v± ex change their signs across the X-point. Note 
that v± ex overshoots v± ix on both two sides of the X-point. 
The outward electron velocity is even faster than the typi- 
cal upstream Alfven velocity. The upstream magnetic field 
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Figure 1. (Color online) Geotail observation from 
1053:00 UT to 1058:30 UT on 15 May 2003. (a) Normal 
magnetic field B z [nT], (b) electric field E x , E y [mV/m], 
and (c,d) plasma perpendicular velocities [km/s] in the 
Xsc and Ysc directions. A potential X-point (red dotted 
line) and a characteristic interval from 1055:07 to 1056:44 
UT (shadow) are overplotted. 



(10 nT) and density (0.01cm" 3 ) give 2200 km/s. Such bi- 
directional, super-Alfvenic electron jets are consistent with 
previous simulations [Karimabadi et al., 2007; Shay et al., 
2007]. 

The shadows indicate a characteristic interval from 
1055:07 to 1056:44 UT, surrounded by sudden spikes in B z 
(Fig. la). This interval features super-Alfvenic electron jets 
(Fig. lc), an intense duskward electron flow (Fig. Id), and a 
low plasma density (not shown). Nagai et al. [2011] referred 
to this interval as the "ion-electron decoupling time." 

Next we introduce several measures to discuss the recon- 
nection region from Geotail data. We consider the energy 
transfer in the rest frame of the electron's bulk flow. This 
measure, the electron- frame dissipation measure D e , is for- 
mally defined in Eq. (7) in Zenitani et al. [2011a]. In a 
nonrelativistic neutral plasma, it can be reduced to 

D e = jE' =j.E-W, (1) 

where E' = E + v e x B is the nonideal electric field and 

W=axB). Vi = (jxB).v e . (2) 

is the work rate per unit volume by the Lorentz force to 
the ion fluid or the electron fluid. Since this Lorentz work 
W (/sec /unit volume) stands for the energy transfer in 
the ideal MHD, Equation 1 reads the total energy transfer 
minus the ideal transfer, i.e., the nonideal energy transfer 
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Figure 2. (Color online) The approximate dissipation 
measure (gray histogram) and the work rate by the 
Lorentz force W (dashed histogram). The X-point (red 
dotted line) and the interval (shadow) are the same as 
those in Figure 1. 



from fields to plasmas [Birn and Hesse, 2005; Zenitani et 
al., 2011a]. We further assume E' z = to obtain the follow- 
ing approximate form, 



Dl = j x E' x + jyE' y . 



(3) 



The two measures and W are shown in Figure 2. The 
current density is directly calculated from the plasma mo- 
ments, j = eni(vi — v e ). The dissipation D* (gray his- 
togram) has a peak near the X-point. Its peak rates are 39 
and 45 pWm" 3 during 1055:32 - 1055:44 UT and 1055:44 - 
1055:56 UT (two 12 s intervals), respectively. The Lorentz 
work W (dashed histogram) has a peak during 1054:31- 
1054:43 UT. 

3. PIC simulation 

We validate our observations with a 2D PIC simula- 
tion [Zenitani et al., 2011b]. We employ a Harris-like ini- 
tial model, B(z) = Bo tanh(jz/L)x and n(z) = rio[0.2 + 
cosh _2 (z/L)], where -Bo and no are the reference magnetic 
field and density. The current sheet thickness is set to 
L = 0.5di, where di is the ion inertial length. Other 
simulation parameters are m,/m E = 100, uj pe /Q ce — 4, 
Ti/T e = 5, and 2400 x 1600 cells with 2.2x10^ particles. 
The computational domain is x, z € [0,76.8] x [—19.2, 19.2] 
in units of di with periodic (x) and reflecting (z) bound- 
aries. Velocities are normalized by the typical Alfven veloc- 
ity c A = B /{^on mi) 1/2 . 

Figure 3 shows snapshots in the well-developed stage. In 
this case, judging from the sign of D e , the dissipation region 
is located around 34.7 < x < 41.5 and -0.38 < z < +0.38. 
Its thickness is limited by that of an electron current layer. 
The layer thickness is comparable with the local electron 
inertial length, d l ° c ~ 0.31. There are also weak dissipative 
regions at x — 28 and x = 48, where the magnetic field lines 
suddenly flip. Surprisingly, both D e and D* give almost 
identidal pictures, as shown in Figures 3a and 3b. There 
are minor differences near the vertical dissipative regions 
but it is hard to recognize them in the figures. 

Figure 3c presents plasma outflow velocities along the 
outflow line (z — 0). The electron jet substantially over- 
runs the ion flow [Karimabadi et al., 2007; Shay et al., 2007] 
until it reaches a shock-like transition region at x w 28 and 
x ~ 48. The electrons are magnetized further downstream. 
Note that Vi X and v ex are still different in the downstream, 
because the ions remain unmagnetized. 

Figures 3d shows the dissipation measure D e and the 
Lorentz work W. The D e measure is almost identical to 
along the outflow line, where j z E' z m 0. Outside the 
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central dissipation region, D e exhibits a weak undershoot in 
the downstream, indicating energy transfer from plasmas to 
the electromagnetic fields in the comoving frame of plasmas. 
This is a signature of the fast electron jet region [Zenitani 
et al., 2011b]. The Lorentz work W is usually positive over 
the reconnection region, because the Lorentz force contin- 
uously drives plasmas in the outflow directions. Further 
downstream, it has a local peak at x«16.5 (outside the do- 
main in Fig. 3), where the reconnected flux (B z ) hits the 
initial current sheet. 

4. Discussion 

Let us examine the validity of the approximate measure 
(Eq. 3). Since EFD does not measure E z along the spin 
axis, we have to reconstruct E z in Equation 1. We consider 
the following three ways. (1) A popular way is to assume 
E ■ B = 0. However, since there are strong parallel electric 
field near the dissipation region, Em 7^ [Pritchett, 2001; 
Wygant et al., 2005], and since B z is usually weak, this 
does not work. (2) The second is to assume E z — 0. This 
gives a different picture that emphasizes separatrices. Phys- 
ically, this assumption drops the energy transfer j z E z but 
leave the ideal part j z (y e x B) z in Equation 1. In other 
words, the ideal part contaminates the nonideal dissipation 
measure D e . Thus we should not assume E z = 0. (3) The 
third is to set E' z — 0. This is our choice. Despite strong 
Hall electric field E z near the dissipation region [Shay et al., 
1998; Wygant et al., 2005], the E' z » condition is fairly 
satisfied outside the electron current layer, and so j z E' z is a 
minor contributor to D e along the inflow line (see Fig. 5a in 
Zenitani et al. [2011b]). The ideal condition is also violated 
E' z 7^ along the separatrices but this does not involve sig- 
nificant energy transfer. In general, j z E' z is negligible due 
to weak j z . Thus fl* approximates D e quite well. Note 
that this choice is specific to a magnetotail configuration 
without a guide field. In other configurations j z E' z could be 
important and Equation 3 may not be useful. 

In Figure 2, the D* measure has a peak near the X-point 
during 1055:44 - 1056:08 UT. This is evident for the dis- 
sipation region. Let us evaluate the dissipation rate. The 
reconnection electric field E y is fairly uniform over the re- 
connection region, except for the magnetic pile-up region 
near the outflow jet front. From PIC simulations, we em- 
pirically know that the reconnection rate remains constant, 
E y ~ 0.1cA,upB up . Here, B up and ca,u P are the magnetic 
field and the Alfven velocity in an upstream region, typi- 
cally a few di upstream from the X-point. It is empirically 
known that the half width of the electron current layer is 
approximated by the local electron inertial length d° c up to 
the case of a realistic mass ratio [Pritchett, 2010]. The cur- 
rent density inside the dissipation region is j ~ Bd/ (nod}° c ), 
where Bd is the magnetic field at the upstream edge of the 
dissipation region. After the initial current sheet plasma 
moves out, the plasma density is roughly uniform over the 
reconnection region. Assuming a uniform density, one can 
estimate typical dissipation at the X-point, 

d: - " (^)( I "~ S -) 

From simulations, we know that the upstream magnetic field 
is typically a half of the asymptotic lobe field : B up ~ \Bq. 
In this event, the lobe field is estimated to Bo = 20 nT, while 
B x occasionally hits 10 nT near the separatricies. Thus, it 
would be reasonable to assume B up = 10 nT. Meanwhile, 
it is very difficult to estimate Bd, because the electron-scale 
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Figure 3. (Color online) (a) The dissipation measure 
D e (Eq. 1) in unit of JqcaBo with magnetic field lines 
superposed, (b) the approximate dissipation measure D* e 
(Eq. 3) in the same unit, (c) Vi X ,v ex at z = 0, and (d) 
the dissipation measure D e and the Lorentz work W at 
z = 0. 



dissipation region is embedded in a thick reconnection layer 
of the order of the ion meandering width. From simulation 
we obtain (Bd/B up ) ~ 0.5. However, this ratio could be 
smaller in a real system, because the electron skin depth 
with rrn/m e = 1836 is 4.3 times smaller than one with 
rm/me = 100. Using B up = 10 nT and (B d /B up ) < 0.5, 
we obtain D* < 163 pWm -3 . In addition, our data are 
averaged over a time interval of 12 s. Geotail resolves the 
spatial structure rather coarsely. Using spatially averaged 
quantities, their product D* could be underestimated by a 
factor of two or three. Considering this, the observed rates, 
39 and 45 pWm -3 (1055:44 - 1056:08 UT), are reasonable. 

It is difficult to evaluate unambiguously the size of the 
dissipation region in our single-satellite observation. Shay 
et al. [2007] predicted that the half length of the E' y > 
region is ~0.6di. In many rrii/m e = 25 runs [A. Klimas, 
private communication] and in our m,i/m e — 100 run, the 
dissipation region is twice longer than the E' y > region. 
Thus the full length of the dissipation region will be 2Adi. 
On the other hand, Zenitani et al. [2011b] reported that the 
aspect ratio of the dissipation region is ~ 0.1. Given that 
the half width of the dissipation region is ~ d 1 ° c , the full 
length will be ~ 20d 1 ° c ~ 0.47d'° c - In these cases, the full 
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extent of the dissipation region will be 

2Ad * ~ 1726 (cuSfO km ' (5) 

0.47d oc ~ 1061 ( — -p ? ) km. (6) 

V 0.01cm- 3 / v ; 

Examining several clues such as the upstream distribution 
functions, Nagai et al. [2011] estimated that the entire struc- 
ture retreats at ~100 km/s. The above two estimates are 
comparable with the spatial scale of one or two sampling 
intervals (12-24 s), 1200-2400 km. Geotail detected the dis- 
sipation region by two or three sampling intervals (Fig. 2), 
2400-3600 km or 3-5cfc or l-2d' oc . The length of the dissi- 
pation region is on the same order. 

The characteristic interval from 1055:07 to 1056:44 UT 
(the shadow regions in Figs. 1 and 2) would correspond to 
the full extent of the outer electron jet region in the PIC 
simulation (28 < x < 48 in Fig. 3). One can recognize 
the magnetic field fluctuations in Figure la (indicated by 
arrows). It is likely that the fluctuations originate from the 
electron jet front region at x ~ 28 and x « 48 (Fig. 3), 
where the fast electron jet hits the outer plasma outflow. 
In observation, D* is negative before the central dissipation 
region arrives (Fig. 2). This is consistent with the weak 
undershoot of D e in the electron jet region in the PIC simu- 
lation (Fig. 3). We do not recognize a similar undershoot on 
the other side (1055:44 - 1056:44 UT). This is because the 
Geotail was near the separatrix, outside the central electron 
jet region [Nagai et al., 2011]. It is not surprising that v±i 
still differs from v± e outside the interval (Figs, lc and Id). 
Electrons are magnetized there, but ions remain unmagne- 
tized downstream of the electron jet fronts. Reconnection 
outflow recovers the ideal MHD approximation when ions 
are magnetized further downstream, but the transition to 
the MHD region could be ambiguous. 

As seen in Figure 2, the Lorentz work W is usually pos- 
itive over the reconnection region. This signature is con- 
sistent with the PIC simulation (Section 3). One can rec- 
ognize a peak during 1054-1055 UT. We find that B z is 
compressed (Fig. la) in this interval and that the plasma 
density decreased after the W-region arrived (see Fig. lj in 
Nagai et al. [2011]). This IT-region appears to be a mag- 
netic pile-up region behind a tangential discontinuity, where 
the compressed B z pushes the dense plasma sheet outward. 
We do not recognize similar W^-region on the other side of 
1057-1058 UT, because the satellite was not in the central 
neutral plane. In fact, we find similar W-regions associated 
with the pile-up regions in other reconnection events. The 
W measure in reconnection deserves further investigation by 
theories, simulations, and observations. 

To our knowledge, this is the first identification of the 
dissipation region in magnetotail reconnection. Earlier ob- 
servations focused on large-scale signatures of the reconnec- 
tion site, such as the flow reversal or Hall fields. In the 15 
May 2003 event, Geotail fortunately resolved a small-scale 
signatures of bi-directional electron jets [Nagai et al., 2011] 
and then the present work finally detected a compact dis- 
sipation region between the jets. Although our resolution 
is quite limited, this is the best possible measurement with 
the present instruments. In addition, this work is a step 
forward for reconnection theories. The dissipation measure 
theory [Zenitani et al., 2011a] was only tested by 2D PIC 
simulations with artificial mass ratio. Here we showed that 
the theory works in the complex real world. The theory 
appears to be practically useful. 

5. Conclusion 

We analyzed Geotail reconnection event on 15 May 2003. 
We applied the electron-frame dissipation measure D e {D*) 
to the reconnection structure. Using D*, for the first time, 



we identified the dissipation region around the possible X- 
point in nature. The dissipation rate (45pWm -3 ) and the 
length of the dissipation region (l-2d' oc ) are reasonable. We 
introduced the Lorentz work W as well. The two measures 
are useful for better understanding of the reconnection struc- 
ture. 

Our approach will be applicable to next-generation in situ 
observations. NASA's upcoming Magnetospheric Multiscale 
(MMS) mission features high temporal resolution, multi- 
satellite observation, and measurement of all three compo- 
nents of the electric field. We will be able to see clearer 
pictures of the dissipation region. 

Acknowledgments. S.Z. acknowledges S. Imada, M. 
Hoshino, Y. Miyashita, M. N. Nishino, and M. H. Saito for useful 
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